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Abstract. We review the nuclear forces currently in use, i. e., the high-precision 
NN potentials of the 1990’s and the nuclear two- and many-body forces based upon 
chiral effective field theory (EFT). We argue that the EFT approach is superior 
to any of the older schemes. Since accurate chiral forces are available now, the 
stage is set for microscopic nuclear struture to move into a new and exciting era. 


1. Introduction 

This workshop deals essentially with microscopic nuclear structure which has the goal 
to derive the properties of atomic nuclei from the ‘elementary’ forces between nucleons. 
Thus, the input for microscopic nuclear structure calculations are the ‘basic’ nuclear 

forces: two-nucleon forces (2NF), three-nucleon forces (3NF),_ It is the purpose 

of this contribution to review the nuclear forces that are presently in use. We can 
distinguish between two groups of current forces: the so-called high-precision nucleon- 
nucleon (NN) potentials that were developed in the previous decade and NN potentials 
(plus associated 3NFs) based upon chiral perturbation theory. Presently, the majority 
of nuclear structure calculations are still conducted with high-precision potentials, but 
the interest in the chiral forces is increasing. The future belongs to the chiral approach 
for reasons that we will explain in sect. 3. 

Diversity has always been a characteristic feature of the NN potential market. 
Therefore, it is not possible to squeeze all current potentials into the above two 
categories. Examples for such exceptional cases are the Moscow potential [TJ , the 
highly nonlocal potentials developed by Doleschall and coworkers |2j, and the NN 
potentials from inverse scattering by Shirokov et al. J3J. All these potentials have 
recently been applied in nuclear structure yielding remarkable results. Unfortunately, 
because of lack of space, I cannot further elaborate on them. 

2. Historical perspective: The high-precision NN potentials 

In 1993, the Nijmegen group [1 published a phase shift analysis which described the 
NN data below 350 MeV laboratory energy available at the time with the ‘perfect’ 
X 2 /datum of 1.0. This raised the expectation that also NN potentials should reproduce 
the NN data with similar precision. Note that even the best NN models of the 1980’s, 
like the Paris [2| and the Bonn 0 potentials, fit the NN data typically with a y 2 /datum 
« 2 or more. To put microscopic nuclear structure theory to a reliable test, one needs 
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Table 1 . % 2 /datum for the reproduction of the 1992 and 1999 NN databases 
below 350 MeV by the Nijmegen phase shift analysis [J] and two high-precision 
potentials: the CD-Bonn potential TO and the Argonne Vis potential 



CD-Bonn 

Nijmegen 

Argonne 


potential 

PSA 

Vis pot. 

proton- 

proton data 



1992 pp database (1787 data) 

1.00 

1.00 

1.10 

After-1992 pp data (1145 data) 

1.03 

1.24 

1.74 

1999 pp database (2932 data) 

1.01 

1.09 

1.35 

neutron 

-proton data 



1992 np database (2514 data) 

1.03 

0.99 

1.08 

After-1992 np data (544 data) 

0.99 

0.99 

1.02 

1999 np database (3058 data) 

1.02 

0.99 

1.07 

pp and np data 



1992 NN database (4301 data) 

1.02 

0.99 

1.09 

1999 NN database (5990 data) 

1.02 

1.04 

1.21 


a perfect NN potential such that discrepancies in the predictions cannot be blamed 
on a bad fit of the NN data. 

So, the Nijmegen analysis triggered a feverish activity among groups traditionally 
involved in potential construction. The output was a new family/generation of NN 
potentials which eventually were dubbed the high-precision potentials. The research 
groups involved and the names of their new creations are, in chronological order: 

• Nijmegen group J7j: Nijm-I, Nijm-II, and Reid93 potentials. 

• Argonne group |S]: Vis potential. 

• Bonn group unnj: CD-Bonn potential. 

All these potentials have in common that they are charge-dependent which is crucial 
to obtain a good fit of the pp and np data. Moreover, they all need about 40-50 
parameters to reproduce the 1992 Nijmegen NN data base with a x 2 /datum « 1. 
Note, however, that since 1993 the pp database has substantially expanded and for 
the current database the x 2 /datum produced by some of these potentials is not so 
perfect anymore (cf. table 1). 

Concerning the theoretical basis of these potential, one could say that they are 
all—more or less—constructed ‘in the spirit of meson theory’ (e.g., all potentials 
include the one-pion-exchange (OPE) contribution). However, there are considerable 
differences in the details leading to considerable off-shell differences among the 
potentials. 

The CD-Bonn potential uses the full, original, nonlocal Feynman amplitude for 
OPE, while all other potentials apply local approximations. As a consequence, the CD- 
Bonn potential has a weaker tensor force as compared to all other potentials. This is 
reflected in the predicted D-state probabilities of the deuteron, Pd, which is a measure 
of the strength of the nuclear tensor force. While CD-Bonn predicts Pd = 4.85%, the 
other potentials yield Pd = 5.7(1)%. These differences in the strength of the tensor 
force lead to considerable differences in nuclear structure predictions. The CD-Bonn 
potentials predicts 8.00 MeV for the triton binding energy, while the local potentials 
predict only 7.62 MeV. We note that the recent highly nonlocal Doleschall potential [2, 
predicts Pd = 3.6% and reproduces the entire triton binding of 8.48 MeV from the 
2NF alone. 
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The OPE contribution to the nuclear force essentially takes care of the long- 
range interaction and the tensor force. In addition to this, all models must describe 
the intermediate and short range interaction, for which very different approaches are 
taken. The CD-Bonn includes (besides the pion) the vector mesons p(769) and w(783), 
and two scalar-isoscalar bosons, a, using the full, nonlocal Feynman amplitudes for 
their exchanges. Thus, all components of the CD-Bonn are nonlocal and the off-shell 
behavior is the original one as determined from relativistic field theory. 

The models Nijm-I and Nijm-II are based upon the Nijmegen78 potential ED 
which is constructed from approximate one-boson-exchange (OBE) amplitudes. 
Whereas Nijm-II uses the local approximations for all OBE contributions, Nijm-I 
keeps some nonlocal terms in the central force component (but the Nijm-I tensor force 
is local). However, nonlocalities in the central force have only a very moderate impact 
on nuclear structure. If one wants to retain nonlocality, it is more important to keep 
the tensor force nonlocalities. 

The Reid93 {7 and Argonne Vis 0 potentials do not use meson-exchange for 
intermediate and short range; instead, a phenomenological parametrization is chosen. 
The Argonne Vis uses local functions of Woods-Saxon type, while Reid93 applies local 
Yukawa functions of multiples of the pion mass, similar to the original Reid potential of 
1968 JT2J. At very short distances, the potentials are regularized either by exponential 
(Vis, Nijm-I, Nijm-II) or by dipole (Reid93) form factors, which are all local functions. 

Over the past ten years, the family of high-precision potentials has provided a 
useful service to the community. Practitioners in the held of microscopic nuclear 
structure and exact few-body calculations could finally be sure that their predictions, 
good or bad, had nothing to do with a bad fit of the NN data. Moreover, since 
all potentials reproduce the NN data base with the same y 2 /datum « 1 (‘phase- 
equivalent’ potentials), the impact of off-shell differences between potentials on nuclear 
structure predictions could be studied systematically and reliably E3 El ESI EE] 

However, in spite of these great practical achievements, the high-precision 
potentials cannot be the end of the story. From a fundamental point of view, these 
potentials are not satisfactory at all. To achieve the acclaimed accuracy, the potentials 
use about 45 parameters. The underlying one-boson-exchange model has only about 
a dozen parameters ED Additional parameters are introduced, e.g., by using partial 
wave dependent coupling constants, which is hard to justify by the underlying model. 
The motto is simply: 

If you want more more accuracy, 
you have to use less theory. 

This is a slap in the face of theoretical physics. Moreover, fundamental questions are 
unanswered, like: 

• What is the connection to the fundamental theory of strong interactions, QCD? 

• Why is the 2NF so much stronger than the 3NF? 

In the next section, we will show how to overcome these problems. 

3. The EFT approach to nuclear forces 

Quantum chromodynamics (QCD), the generally accepted theory of strong 
interactions, is non-perturbative in the low-energy regime characteristic for nuclear 
physics. For many years, this fact was perceived as the great obstacle for a derivation 
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of the nuclear force from QCD—impossible to overcome except by lattice QCD. The 
effective field theory (EFT) concept has shown the way out of this dilemma. Notice 
that the QCD Lagrangian for massless up and down quarks is chirally symmetric, i. e., 
it is invariant under global flavor SU(2)l x SU(2)r equivalent to SU(2)y x SU(2)a 
(vector and axial vector) transformations. The axial symmetry is spontaneously 
broken as evidenced in the absence of parity doublets in the low-mass hadron spectrum. 
This implies the existence of three massless Goldstone bosons which are identified with 
the three pions (tt^,^ 0 ). The non-zero, but small, pion mass is a consequence of the 
fact that the up and down quark masses are not exactly zero (some small, but explicit 
symmetry breaking). Thus, we arrive at a low-energy scenario that consists of pions 
and nucleons interacting via a force governed by spontaneously broken approximate 
chiral symmetry. 

The effective Lagrangian describing this scenario is given by an infinite series of 
terms with increasing number of derivatives and/or nucleon fields, with the dependence 
of each term on the pion field prescribed by the rules of broken chiral symmetry. 
Applying this Lagrangian to NN scattering generates an unlimited number of Feynman 
diagrams. However, Weinberg showed m that a systematic expansion exists in terms 
of (Q/A x ) u , where Q denotes a momentum or pion mass, A x « 1 GeV is the chiral 
symmetry breaking scale, and is > 0 (cf. figure 1). This has become known as chiral 
perturbation theory (yPT). For a given order is, the number of contributing terms is 
finite and calculable; these terms are uniquely defined and the prediction at each order 
is model-independent. By going to higher orders, the amplitude can be calculated to 
any desired accuracy. Thus, the motto is now: 

If you want more more accuracy, 

you have to use more theory (more orders). 

This sounds more like what we expect from theoretical physics. Moreover, we have 
connected with QCD. 

Following the first initiative by Weinberg GB], pioneering work was performed by 
Ordonez, Ray, and van Kolck mm who constructed a NN potential in coordinate 
space based upon yPT at next-to-next-to-leading order (NNLO; is = 3). The results 
were encouraging and many researchers became attracted to the new field. Kaiser, 
Brockmann, and Weise m presented the first model-independent prediction for the 
NN amplitudes of peripheral partial waves at NNLO. Epelbaum et al. J23 developed 
the first momentum-space NN potential at NNLO, and Entem and Machleidt m 
presented the first potential at N 3 LO. 

In xPT, the NN amplitude is uniquely determined by two classes of contributions: 
contact terms and pion-exchange diagrams. There are two contacts of order Q° 
[O(Q 0 )] represented by the four-nucleon graph with a small-dot vertex shown in the 
first row of figure 1. The corresponding graph in the second row, four nucleon legs and 
a solid square, represent the seven contact terms of 0(Q 2 ). Finally, at 0(Q 4 ), we have 
15 contact contributions represented by a four-nucleon graph with a solid diamond. 

Now, turning to the pion contributions: At leading order [LO, O(Q 0 ), is = 0], 
there is only the wellknown static one-pion exchange (OPE), second diagram in the 
first row of figure 1. Two-pion exchange (TPE) starts at next-to-leading order (NLO, 
is = 2) and all diagrams of this leading-order two-pion exchange are shown. Further 
TPE contributions occur in any higher order. Of this sub-leading TPE, we show only 
two representative diagrams at NNLO and three diagrams at N 3 LO. While TPE at 
NNLO was known for a while (HI Et H2], TPE at N 3 LO has been calculated only 
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2N Force 3N Force 4N Force 




Figure 1. Hierarchy of nuclear forces in yFT. Solid lines represent nucleons and 
dashed lines pions. Further explanations are given in the text. 


recently by Kaiser [20. All 27r exchange diagrams/contributions up to N 3 LO are 
summarized in a pedagogical and systematic fashion in Ref. 2 where the model- 
independent results for NN scattering in peripheral partial waves are also shown. 

Finally, there is also three-pion exchange, which shows up for the first time at 
N 3 LO (two loops; one representative 37 t diagram is included in figure 1). In Ref. 20 j 
it was demonstrated that the 37r contribution at this order is negligible. 

One important advantage of \PT is that it makes specific predictions also for 
many-body forces. For a given order of yPT, 2NF, 3NF, ... are generated on the 
same footing (cf. figure 1). At LO, there are no 3NF, and at next-to-leading order 
(NLO), all 3NF terms cancel |18l I27|. However, at NNLO and higher orders, well- 
defined, nonvanishing 3NF occur 123120 Since 3NF show up for the first time at 
NNLO, they are weak. Four-nucleon forces (4NF) occur first at N 3 LO and, therefore, 
they are even weaker. 
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Table 2. Number of parameters used for the np potential in various approaches 
discussed in the text. 



Nijmegen 

PWA93 

CD-Bonn 

“high 

precision” 

NLO 

Q 2 

(NNLO) 

N 3 LO 

Q 4 


3 

4 

2 

4 

3 Si 

3 

4 

2 

4 

3 Si- 3 £>i 

2 

2 

1 

3 

L Pi 

3 

3 

1 

2 

3 Po 

3 

2 

1 

2 

3 Pi 

2 

2 

1 

2 

3 P 2 

3 

3 

1 

2 

■ i P2~ 3 F2 

2 

1 

0 

1 

l d 2 

2 

3 

0 

1 

3 Di 

2 

1 

0 

1 

3 d 2 

2 

2 

0 

1 

3 d 3 

1 

2 

0 

1 

6 d 3 - 6 g 3 

1 

0 

0 

0 

l f 3 

1 

1 

0 

0 

3 F2 

1 

2 

0 

0 

3 f 3 

1 

2 

0 

0 

3 f 4 

2 

1 

0 

0 

i F4- i H 4 

0 

0 

0 

0 

"G4 

1 

0 

0 

0 

3 g 3 

0 

1 

0 

0 

3 g 4 

0 

1 

0 

0 

3 g 5 

0 

1 

0 

0 

Total 

35 

38 

9 

24 


4. Chiral NN potentials 

The two-nucleon system is non-perturbative as evidenced by the presence of shallow 
bound states and large scattering lengths. Weinberg CHI showed that the strong 
enhancement of the scattering amplitude arises from purely nucleonic intermediate 
states. He therefore suggested to use perturbation theory to calculate the NN 
potential and to apply this potential in a scattering equation (Lippmann-Schwinger 
or Schrodinger equation) to obtain the NN amplitude. We follow this philosophy. 

Chiral perturbation theory is a low-momentum expansion. It is valid only for 
momenta Q <C A x « 1 GeV. Therefore, when a potential is constructed, all expressions 
(contacts and irreducible pion exchanges) are multiplied with a regulator function, 



where p and p' denote, respectively, the magnitudes of the initial and final nucleon 
momenta in the center-of-mass frame; and A <C A x . The exponent 2 n is to be chosen 
such that the regulator generates powers which are beyond the order at which the 
calculation is conducted. 

NN potentials based upon yPT at NNLO 22' 2U1 are poor in quantitative terms; 
they reproduce the NN data below 290 MeV lab. energy with a y 2 /datum of more 
than 20 (cf. tables 3 and 4, below). As shown first by Entern and Machleidt m , one 
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Table 3. x 2 /datum for the reproduction of the 1999 np database below 290 
MeV by various np potentials. (A = 500 MeV in all chiral potentials.) 


Bin (MeV) 

# of data 

N 3 LO“ 

NNLO 6 

NLO 6 

AV18 C 

0-100 

1058 

1.06 

1.71 

5.20 

0.95 

100-190 

501 

1.08 

12.9 

49.3 

1.10 

190-290 

843 

1.15 

19.2 

68.3 

1.11 

0-290 

2402 

1.10 

10.1 

36.2 

1.04 


“Reference na b Reference 122. “Reference 0. 


has to go to order N 3 LO to obtain a NN potential of acceptable accuracy. Therefore, 
we will discuss now specifically the NN potential at N 3 LO. 

At N 3 LO, there are 24 contact terms (24 parameters) which contribute to the 
partial waves with L < 2. In table 2, column ‘N 3 LO/<5 4 ’, we show how these 
terms/parameters are distributed over the various partial waves. For comparison, 
we also show the number of parameters used in the Nijmegen partial wave analysis 
(PWA93) 0 and in the high-precision CD-Bonn potential m- The table reveals that, 
for S and P waves, the number of parameters used in high-precision phenomenology 
and in EFT at N 3 LO are about the same. Thus, the EFT approach provides 
retroactively a justification for what the phenomenologist of the 1990’s were doing. 
At NLO and NNLO [3171 . the number of parameters is substantially smaller than for 
PWA93 and CD-Bonn, which explains why this order is insufficient for a quantitative 
potential. 

For an accurate fit of the low-energy pp and np data, charge-dependence is 
important. Charge-dependence up to next-to-leading order of the isospin-violation 
scheme (NL0, in the notation of Ref. ED) includes: the pion mass difference in OPE 
and the Coulomb potential in pp scattering, which takes care of the L0 contributions. 
At order NLO we have pion mass difference in the NLO part of TPE, 7T7 exchange E2j, 
and two charge-dependent contact interactions of order Q° which make possible an 
accurate fit of the three different 1 So scattering lengths, a pp , a nn , and a np . 

In the optimization procedure, we fit first phase shifts, and then we refine the fit 
by minimizing the y 2 obtained from a direct comparison with the data. The y 2 /datum 
for the fit of the np data below 290 MeV is shown in table 3, and the corresponding 
one for pp is given in table 4. The y 2 tables show the quantitative improvement of the 
NN interaction order by order in a dramatic way. Even though there is considerable 
improvement when going from NLO to NNLO, it is clearly seen that N 3 LO is needed 


Table 4. x 2 /datum for the reproduction of the 1999 pp database below 290 
MeV by various pp potentials. (A = 500 MeV in all chiral potentials.) 


Bin (MeV) 

# of data 

N 3 LO“ 

NNLO 6 

NLO 6 

AV18 C 

0-100 

795 

1.05 

6.66 

57.8 

0.96 

100-190 

411 

1.50 

28.3 

62.0 

1.31 

190-290 

851 

1.93 

66.8 

111.6 

1.82 

0-290 

2057 

1.50 

35.4 

80.1 

1.38 


“Reference ESI. b See footnote 1331 . Reference 0. 










EFT of two- and many-body forces 


3 75 F 

S’ * 
5. 50 

1 25 
to 

8> o 
£ 

-25 


1 Sr 







0 100 200 300 

Lab. Energy (MeV) 


0 100 200 300 

Lab. Energy (MeV) 


0% 3 Ph 


W 


10 

3S -2“ 
(0 

£ -30 




^.l 

l , | i __js 1 

0 100 200 300 

Lab. Energy (MeV) 



0 100 200 300 

Lab. Energy (MeV) 


0 100 200 300 

Lab. Energy (MeV) 



0 100 200 300 

Lab. Energy (MeV) 




Lab. Energy (MeV) 


Lab. Energy (MeV) 




Figure 2. np phase parameters below 300 MeV lab. energy for partial waves with 
J < 2. The thick solid (dashed) line is the result by Entem and Machleidt at 
N 3 LO using A = 500 MeV (A = 600 MeV). The thin dotted and dash-dotted lines 
are the phase shifts at NLO and NNLO, respectively, as obtained by Epelbaum 
et al. I M using A = 500 MeV. The solid dots show the Nijmegen multienergy np 
phase shift analysis [||, and the open circles are the GWU/VPI single-energy np 
analysis SM99 1351 . 


to achieve an accuracy comparable to the phenomenological high-precision Argonne 
Vis potential 0- Note that proton-proton data have, in general, smaller errors than 
np data which explains why the pp x 2 are always larger. 

The phase shifts for np scattering below 300 MeV lab. energy are displayed in 
figures 2. What the \ 2 tables revealed, can be seen graphically in this figure. The 
2 P 2 phase shifts are a particularly good example: NLO (dotted line) is clearly poor. 
NNLO (dash-dotted line) brings improvement and describes the data up to about 100 
MeV. The difference between the NLO and NNLO curves is representative for the 
uncertainty at NLO. and, similarly, the difference between NNLO and N 3 LO reflects 
the uncertainty at NNLO. Obviously, at N 3 LO (A = 500 MeV, thick solid line) we 
have a good description up to 300 MeV. An idea of the uncertainty at N 3 LO can be 
obtained by varying the cutoff parameter A. The thick dashed line is N 3 LO using 
A = 600 MeV. In most cases, the latter two curves are not distinguishable on the scale 
of the figures. Noticeable differences occur only in 1 D 2 , 3 -F2, an d e 2 above 200 MeV. 

5. Chiral three-nucleon forces 

As noted before, the first non-vanishing 3NF terms occur at NNLO and are shown in 
figure 1 (row ‘Q 3 /N 2 LO’, column ‘3N Force’). There are three diagrams: the TPE, 
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OPE, and 3N-contact interactions 

2 

a3NF 


r/di\r _ 

'TPE — 


9A 
2 fn 


. The TPE 3N-potential is given by 

(<7z ■ $)(vj ■ qj) 


( q 2 


ml)(q? +ml) 


F a0 a 0 

r ijk i j 


( 2 ) 


with ifi = pf — pi , where pi and pf are the initial and final momenta of nucleon i, 
respectively, and 


Ft = 600 


4cim‘ 2 c 3 


f 2 

J 7T 


~f2 ® 
J 7T 


£i \ ' c «/37 _7 


■'7T ~ 


fit cr k ■ [qi x qj] .(3) 


The vertex involved in this 3NF term is the two-derivative ttttNN vertex (large solid 
dot in figure 1) which we encountered already in the TPE contribution to the 2N 
potential at NNLO. Thus, there are no new parameters and the contribution is fixed 
by the LECs used in NN. The OPE contribution is 

9a \ ' &j ■ qj 
8 ^ qf + ml 

i^j^k n 


y3NF _ n 
'OPE — u 


■ Tj){Si ■ qj) 


(4) 


and, finally, the 3N contact term reads 

Kf F = E\Y, T rr k . (5) 

The last two 3NF terms involve two new vertices (that do not occur in the 2N 
problem), namely, the nNNNN vertex with parameter D and a 6iV vertex with 
parameters E, To pin them down, one needs two 3N observables. In reference j28l . 
the triton binding energy and the nd doublet scattering length 2 a n d were used. Once 
D and E are fixed, the results for other 3N, 4N, ... observables are predictions. In 
reference ESI, encouraging results were reported for 6 Li. Concerning the famous l A y 
puzzle’, the above 3NF terms yield some improvement of the predcited nd A yi however, 
the problem is not resolved. 

One should note that there are additional 3NF terms at NNLO due to relativistic 
corrections (1 /Mjv corrections) that have not yet been included in any calculation. 
However, there are all reasons to believe that these contributions will be very small, 
probably negligible. It is more likely that the problem with the chiral 3NF is analogous 
to the one with the chiral 2NF: namely, NNLO is insufficient and for sufficient accuracy 
one has to proceed to N 3 LO. Two 3NF terms at N 3 LO are indicated in figure 1. The 
N 3 LO 3NF, which will probably not depend on any new parameters, is presently under 
development. 


6. Conclusions 

The EFT approach to nuclear forces is superior to all earlier nuclear force models. It 
represents a scheme that has an intimate relationship with QCD and allows to calculate 
nuclear forces to any desired accuracy. Moreover, nuclear two- and many-body forces 
are generated on the same footing. 

At N 3 LO 22, the accuracy is achieved that is necessary and sufficient for reliable 
nuclear structure calculations. First calculations applying the N 3 LO potential have 

produced promising results (S3 ESI E2 EH OH 113 El 

The 3NF at NNLO is known J2BJ and has had first successful applications 135) . 
The 3NF and 4NF at N 3 LO is presently under construction. 

In summary, the stage is set for many years of exciting nuclear structure research 
that is more consistent than anything we had before. 
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